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The first observations of annihilation processes of the type e + e~ — > J/if)T] c and J /if) (cc) non -res 
were recently reported by Belle [lj]. The measured cross-sections for both processes are 
an order-of-magnitude larger than theoretical predictions based on non-relativistic QCD 
(NRQCD) 0, [|. In an attempt to explain at least part of this discrepancy, the authors 
of Ref. |4j suggest that processes proceeding via two virtual photons may be important. 
In particular, if the two-photon-mediated process e + e~ — > J /if) J /if) has a significant cross- 
section, the observed e + e" — > J/if>r] c signal, which is inferred from the r] c peak in the 
recoil mass spectrum for the reconstructed J /if) in inclusive e + e~ — > J/ipX events, might 
also include double J /if) events and, thus, produce an inflated cross-section measurement. 
e + e~ annihilation to J /if) J /if) via a single virtual photon is forbidden by charge conjugation 
symmetry and was ignored in our published analysis. Here, using a data sample of 101.8 
fb _1 collected by the Belle detector and the analysis procedure described in Ref. jlj, we 
evaluate this possibility. 

Since the r] c and J /if) have similar masses (Mjm — M Vc ~ 116 MeV/c 2 ), it is important to 
check for any momentum scale bias that may shift the recoil mass values. We use e + e _ — > 
if)(2S)j, if)(2S) — > J/if)7T + 7T~ events to calibrate and verify the recoil mass scale. We find 
that any shift in the recoil mass is less than 3 MeV/c 2 . 

The spectrum of recoil masses against the J /if) is presented in Fig. H a clear peak is 
observed around the r] c nominal mass, and a smaller peak is seen around the Xco nominal 
mass; the large peak at ~ 3.63GeV/c 2 is interpreted as the r) c (2S). We performed a fit 
to this spectrum that includes all of the known narrow charmonium states. In this fit, the 
mass positions for the T] c , Xco and i] c (2S) are treated as free parameters; those for the J /if), 
Xci, Xc2 and if>(2S) are fixed at their nominal values. The expected line-shapes for these 
peaks are determined from a Monte Carlo simulation as described in our previous paper [IJ , 
the background is parametrized by a second order polynomial function, and only the region 
below the open charm threshold (M recoi i < 3.7 GeV/c 2 ) is included in the fit. The fit results, 
listed in Table HI give negative yields for the J /if), Xci, Xc2 and if)(2S); the solid line in 
Fig. is the result of a fit with all these contributions fixed at zero. The dotted line in 
the figure corresponds to the case where the contributions of the J/if>, Xci, Xc2 and if)(2S) 
are set at their 90% confidence level upper limit values. The dashed line is the background 
function. To set a conservative upper limit for e + e~ — > J /if) J /if), we use assumptions for the 
production and helicity angle distributions that correspond to the lowest detection efficiency. 

In summary, using a larger data set we confirm our published observation of e + e~ — > 
J/if>r) c and find no evidence for the process e + e~ — > J/if)J/if). We set an upper limit for 
a{e+e- -> J /if) J /if)) x B{J/ip -► > 2 charged) of less than 0.008 pb at the 90% CL. 

Although the limit presented here is not inconsistent with the prediction for the J /if) J /if) 
production rate given in Ref. 0], the suggestion that a significant fraction of the inferred 
J/ipr) c signal is actually J /if) J /if) is ruled out. Therefore, the discrepancy between the Belle 
result and the NRQCD prediction remains. 
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FIG. 1: Distribution of masses recoiling against the reconstructed J/ip in inclusive e + e 
events. The curves are explained in the text. 



TABLE I: Summary of the signal yields, charmonium masses and significances for e + e 
J/tp (cc) res . 



(cc) res state 



N 



M [GeV/c 2 



a 



Vc 

J 111) 
XcO 

Xcl + Xc2 

Vc(2S) 
if>(2S) 



175 ± 23 
-9 ±17 
61 ± 21 
-15 ± 19 
107 ± 24 
-38 ±21 



2.972 ± 0.007 

fixed 
3.409 ± 0.010 

fixed 
3.630 ± 0.008 

fixed 
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